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We explore Loschmidt echo in two regimes of one-dimensional (ID) interacting Bose gases: the 
strongly interacting Tonks- Girardeau (TG) regime, and the weakly-interacting mean-field regime. 
We find that the Loschmidt echo of a TG gas decays as a Gaussian when small perturbations 
are added to the Hamiltonian (the exponent is proportional to the number of particles and the 
magnitude of a small perturbation squared). In the mean-field regime the Loschmidt echo decays 
faster for larger interparticle interactions (nonlinearity) , and it shows richer behavior than the TG 
Loschmidt echo dynamics, with oscillations superimposed on the overall decay. 
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I. INTRODUCTION 

The understanding of why an isolated (interacting many-body) system, which is initially say far from equilibrium, 
in many cases macroscopically undergoes irreversible evolution towards an equilibrium state, despite the fact that the 
microscopic laws are reversible, has intrigued scientists ever since the first disputes between Boltzmann and Loschmidt 
on this topic P, @- In principle, if the time was reversed at a given instance, the system would evolve back into the 
initial state. However, such a reversal is for all practical purposes impossible due to high sensitivity to small errors 
and interaction of the system with the environment. The quantity that measures sensitivity of quantum motion 
to perturbations is called Loschmidt echo or fidelity (f° r a review see e.g. 0). Fidelity tells us what is the 

probability that system will end up in the initial state after forward evolution for time t, followed by the slightly 
imperfect time reversed evolution for the same time t. In quantum mechanics time evolution from an initial state 
■00 is given by the unitary operator U(t) = exp(—iHt/h) via ip(t) = exp(—iHt/h)tpo, and the echo dynamics can be 
formally written as 

F(t) = \(i> Q \eMiH>t/h)eM-iHt/h)\^ Q )\ 2 . (1) 

Here H is the Hamiltonian of the unperturbed system, and H' = H + V E is the slightly perturbed Hamiltonian ||. 
One can think about this quantity as measuring the stability of quantum motion Q, i.e., it tells us the overlap of the 
two states ip(t) and ip'(t), the former is evolved forward in time by H, and the latter by H': 

F(t) = \{yj'(t)\m)\ 2 - (2) 

In quantum systems, depending on the strength of the perturbation and the properties of the nonperturbed Hamilto- 
nian, three different decay regimes of Loschmidt echo are usually identified: the Gaussian perturbative regime 0, 0] , 
the exponential Fermi golden rule regime 043] > and the Lyapunov regime Q ■ 

Motivated by the recent progress in experiments and theory on ultracold atomic gases Q, where the influence 
of the environment can be made very small, and the strength of the atom-atom interactions can be tuned @, we 
^ | are motivated to investigate Loschmidt echo dynamics in those systems. In particular, we focus on one-dimensional 
(ID) Bose gases which were experimentally realized fioj even in the strongly correlated regime of Tonks- Girardeau 
(TG) bosons [n], [l2[ , both in [ll| and out of equilibrium |l2| . The realization of the TG gas in atomic waveguides 
was proposed by Olshanii [l~3|. The ID atomic gases can be described with the Lieb-Liniger model [3, which for 
weak interactions is well-described by the Nonlinear Schrodinger equation (Gross-Pitaevskii theory) [9[|, whereas for 
sufficiently strong interactions one enters the Tonks-Girardeau regime, where exact solutions can be found via Fermi- 
Bose mapping [15(. This method was used to study out-of-equilibrium dynamics in the strongly correlated regime 
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(e.g., see [l6l - |2T| ]). The Loschmidt echo was within the mean- field Gross-Pitaevskii theory addressed in Refs. [22I l23j . 
We would also like to point out at a study of orthogonality catastrophe (and the relation to Loschmidt echo) in an 
ultracold Fermi gas coupled to a single cubit [24j . 

Here we demonstrate, with exact numerical calculation, that for small random stationary perturbation, the 
Loschmidt echo for a TG gas decays as a Gaussian with decay constant proportional to the number of particles 
and the square of the amplitude of the perturbation. We analytically derive the Gaussian behavior of TG fidelity 
within approximation presented by Peres Q. In the mean- field regime the Loschmidt echo decays faster for larger 
interparticlc interactions (nonlincarity) , and it shows richer behavior than the TG Loschmidt echo dynamics, with 
oscillations superimposed on the overall decay. 



II. THE PHYSICAL SYSTEM AND THE CORRESPONDING MODEL 



Consider a gas of N identical bosons in a one-dimensional (ID) space, which interact via pointlike interactions, 
described by the Hamiltonian 



H 



N 

E 



2m dX} 



U(Xi) 



9lD 

l<i<j<N 



(3) 



Such a system can be realized with ultracold bosonic atoms trapped in effectively ID atomic waveguides [10H12J, 
where U(X) is the axial trapping potential, and gio = 2h 2 a3D[ma 2 L (l — Ca^o / V^a^)] -1 is the effective ID coupling 
strength; 0323 stands for the three-dimensional s-wave scattering length, a± = \fhjmuj\_ is the transverse width of 
the trap, and C = 1.4603 [l3j . By varying ui± the system can be tuned from the mean field regime described by 
the Gross-Pitaevskii equation, up to the strongly-interacting Tonks-Girardeau regime (gm — > 00). In equilibrium, 
different regimes of these ID gases are usually characterized by a dimcnsionlcss parameter 7 = mgn>/h , nic, where 
n\D stands for the linear atomic density. For 7 <C 1 the gas is in mean field regime and for 7 3> 1 it is in the strongly 
interacting regime (we consider repulsive interactions 7 > 0) [llHlij . One can tunc 7 by say changing the transverse 
confinement frequency. In our calculations, we use Hamiltonian ([3]) in its dimcnsionless form 



H 



N 

E 



o_ 2 

dx 



N 



2c S(xi — Xj). 



(4) 



where x = X/Xq [Xq is the spatial scale which we choose to be l^m). Here we consider 87 Rb atoms with the 3D 
scattering length a^o = 5.3 11m [TJ [l2j]. Energy is in units of Eq = h 2 /2mX^ = 3.82 • 10~ 32 J, and time is in units 
of To = 2mXp/fi = 2.8 ms. The dimensionless axial potential is V{x) = U(X)/Eq, and the interactions strength 
parameter is 2c = gio / 'XqEq. 



III. LOSCHMIDT ECHO OF A TONKS-GIRARDEAU GAS 

For the Tonks-Girardeau gas, the interaction strength is infinite c — > 00, that is, the bosons are "impenetrable" [l5| . 
Consequently, an exact (static and time-dependent) solution of this model can be written via Girardeau's Fcrmi-Bosc 
mapping |15l.ll6j 

ip B (xi,...,x N ,t)= sgn(xi - Xj)ip F (xi,...,x N ,t), (5) 

l<i<j<N 

where ipp denotes a wave function describing N noninteracting spin polarized fermions in the external potential V(x). 
In our simulations we consider up to N = 70 particles. The system is initially (for times t < 0) in the ground state 
of a container like potential, 

V L (x) = V {1 + t&nh[V s (x - L/2)}/2 - tanh[V s (o; + L/2)]/2}, (6) 

where Vb = 500, V s = 4, and L = 15 (corresponding to 15 /im). At t = we suddenly expand the width of the container 
to twice its original width, that is, the potential at times t > is V2l(x). In order to calculate the fidelity F(t), we must 
evolve the TG gas in the new potential Vq,l{x), and in the potential Vl L {x) = V2l(x) + V e (x) starting from identical 
initial states. Here V E (x) is a small noise potential of amplitude e. The Loschmidt echo F(t) — KV'bWIV'b^))! 2 
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is calculated from the knowledge of the TG many-body states ipB(t) and ip' B (t) corresponding to the evolution in 
potentials V2l(x) and V^ L (x), respectively. 

The fermionic wave function ipp can in our case be written as a Slater determinant, ipp[x\, . . . ,xw,t) = 
det m n= i[4>m{xn, i)]/y/N\, where il> m (x, t), m = 1, . . . , N , satisfy the single-particle Schrddinger equation 



.di> m (x,t) 



2 



& 



1pm(x,t), (7) 



and equivalcntly for ij; m (x,t) which evolve in V^(x). The initial conditions arc such that ip m (x, 0) = ip m (x, 0) is the 
m-th single particle eigenstate of the initial container potential Vl(x). The Loschmidt echo for a Tonks-Girardeau 
gas can be written in a form convenient for calculation: 

CTl ?— 1 CT 2 j — 1 

1 N 

- i^EEc-rc-rn^w^wWi 8 

CTl CT2 i—1 

= |detP(t)| 2 , (8) 
where a denotes a permutation in N indices, (— ) a is its signature, and 

Pij(t) = J ^'*{x,t)iP 3 {x,t)dx. (9) 

In writing relation ([5]) we used a definition of the determinant. Since at t = we have Pjy(0) = that motivates us 
to define the fidelity product 

N 

Fp(t) = J[P ii (t)P* i (t). (10) 

i=l 

Thus, in calculation of the fidelity product we assume that all off diagonal elements of the matrix arc zero i.c 
Pijif) = for i ^ j for all times. It can be interpreted as if we evolve the N particles fully independently of each 
other (including statistics) starting from the N initial states ip m (x, 0), m = 1, . . . , N, calculate N different fidelities 
for these states, and multiply them to obtain the product fidelity. The value F(t) is identical for noninteracting 
spinless fermions and interacting TG bosons; note that Eq. (|5J| is identical to the formula used by Goold et al. in 
Ref. [Hi studying the orthogonality catastrophe for ultracold fermions. Thus, Fp and F will distinguish the influence 
of antysymmetrization in the case of noninteracting fermions, or TG interactions and symmetrization in the case 
of bosons, with dynamics which does not include neither statistics nor interactions into account. We would like to 
emphasize that derivation of ([5]) and © does not require that we initiate the dynamics from the ground state of the 
TG gas in the initial trap; we could have chosen any excited TG eigenstate as initial condition as well. 

In order to calculate the fidelity F(t), we must evolve the single particle states ipj(x,t) [?/>-(x, i), respectively], in 
the potential V2l(x) [V2l(x) + V £ (x)], starting from the first N single-particle cigenstates of Vl(x). The evolution is 
performed via standard linear superposition in terms of the eigenstates </> m (x) of the final container potential V2l{x) 
(which are calculated numerically): 

ipj(x,t) = a? m cj> m {x) exp(-iE m t), (11) 

m 

and 

^(x,t) =J2^'m(x)exp(-iE' m t). (12) 
in 

We numerically calculate the coefficients a 3 n = / ^(x)t/jj (x, 0)dx and a? n = J </>* l (x)-!/>j (x, 0)dx for j = 1, . . . , 70, and 
in = 1, ... , 210 which is sufficient for the parameters we used. The noise potential is constructed as follows: x-space is 
numerically simulated by using 2048 equidistant points in the interval x € [—30, 30]. From this array we construct a 
random array V ian d(x) = l-FT -1 [cxp(— fc 4 //^ ut )i ? T[rand(x)]]| of the same length, where rand(x) is a random number 
between and 1 , FT stands for the Fourier transform, K cut is the cut-off wave vector (set to K cut = 53) introduced 
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to make the discrete numerical potential sufficiently "smooth" from point to point. Finally the noise potential is 
obtained via V E (x) = £[V ia , n ^(x) — KandL where e is the amplitude of the perturbation, and Viand is the mean value 
of V raun d(x). Such a potential can be constructed optically for ID Bose gases plj . 

The fidelity depends on the particle number N, the amplitude of noise e, but also on the particular realization of 
V e (x), hence, we calculate all quantities (e.g., the Loschmidt echo) for 50 different realizations of V e (x), and then 
perform the average over the noise ensemble: (F(t)) no i se - 

In Fig.QJa) we show the fidelity (F(t)) no i se as a function of time, for three different numbers of particles, N = 10, 
20 and 50 with e = 0.05. We find that in the TG regime, the Loschmidt echo decays as a Gaussian: (F(t)) no i se — 
exp(— (A(iV, e))t 2 ); solid black lines represent the Gaussian curves fitted to the numerically obtained values. We point 
out that in every single realization of the noise, the fidelity for a TG gas decays as a Gaussian, with small fluctuations 
in the value of the exponent. The error bars in Fig. [TJa) represent the standard deviation of the fidelity at a given time. 
Note that the standard deviation gets smaller with increasing particle number N, which means that for sufficiently 
large N it suffices to calculate fidelity decay for a single realization of the potential to obtain reliable values for X(N, e). 
It s interesting to compare the fidelity with the fidelity-product {Fp{t)) no ise, which can also be fitted well with the 
Gaussian function as illustrated in Fig.QJb). We find that the product (Fp(t)) n oise is systematically below the value 
of the fidelity. 

In Fig. [2] we depict the dependence of the fidelity, that is, of the exponent (X(N, e)), on the number of particles N 
and e. In Fig. (5Ja) we plot (A(7V, e)) je 2 as a function of e for different values of N; evidently we have (X(N, e)) cx e 2 . 
In Fig. [5Jb) we plot (X(N,e))/N as a function of TV for e = 0.05; we clearly see that (X(N,e)) cx TV for sufficiently 
large N (already for N > 20). 
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FIG. 1: (color online) Decay of the Loschmidt echo (fidelity) with time for e = 0.05. (a) The averaged values {F{i)) n0 ise for 
N = 10 (crosses), N — 20 (asterisks), and N = 50 (circles). Solid black lines represent the Gaussian functions fitted to the 
numerically obtained values. Error bars depict the standard deviation for 50 different realizations of the noise potential V E (x). 
(b) The fidelity (F(t)} no ise (circles), fidelity product (Fp(t)) n0 ise (black dot-dashed lines), the values obtained via det(PP T ), 
where Pij = a"aj, exp(icj„t) is obtained via approximation presented by Peres [j| (blue line), and the fidelity obtained via 
trace-log formula Eq. (|14|) (red dotted line). 



In order to understand numerical results of Fig. Q] and Fig. [2] we analytically explore the properties of the fidelity 
for a single realization of V e (x). To this end we use an approximation from Peres [3(, where to first order in e one has 
J dx(j)*4>i f=a Sij, and a 3 m rj a? m . The elements of the matrix P, which yield the fidelity via Eq. ([5]). are then written 
as Pij = ^2 n a"a J n exp(iui n t), where w„ = E' n — E n (<fin\V e \4> n ) ■ In Fig. QJb) we plot the fidelity obtained with this 
approximation (blue line) and the one obtained with the exact numerical evolution (blue circles); the agreement is 
excellent. The diagonal elements |Pii(t)| 2 can be interpreted as single-particle fidelities corresponding to the initial 
states ipi(x, 0). It is straightforward to see that |Pii(t)| 2 = J2 n m \ a n\ 2 \ a m\ 2 cos [( w n — w m )i], however, we note that 
in our simulations only a few terms contribute to the sum above, yielding oscillatory behavior of the single particle 
fidelities with relatively high amplitudes of the oscillation. Now we turn to the TG gas and our observation that the 
decay of fidelity is Gaussian. In order to derive this we use the trace-log formula for the determinants: 

F(t) = exp(Tr(log(PPt))). (13) 

We can approximate PP 1 w 1 + Qi< - Q2^ 2 + 0(t 3 ), where A„ m = w„ - w m , [Qi]ij = i J2k=i E n , m Caid^m^m. 
and [Q2]ij = \ Sfc=i J2 n m a n a n a m a m ^nm- Next we expand the logarithm in trace-log formula which yields 



F(t) = cxp(-TrQ 2 t 2 ), 



(14) 
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FIG. 2: (color online) The Gaussian exponent of the fidelity as a function of e and N. (a) The quantities (A)/e 2 are plotted 
for different particle numbers, and they are ordered just as in the legend (higher lines are for larger values of TV); obviously 
(A) oc e 2 . (b) The quantities {Xp)/N (red circles, upper line) and {\}/N (blue asterisks, lower line), are plotted as a function 
of N for e = 0.05. For larger TV the lines become horizontal indicating that (A) oc TV oc (Ap). 



i.e., a Gaussian function. In our derivation we used TrQi — 0. Red dotted line in Fig. [T^b) shows that Eq. (|14p 
is an excellent approximation for larger N. The dependence of (X(N,e)) on e follows from the fact that A 2 m oc e 2 , 
whereas TrQ 2 oc N (sec Fig. [2J. 

In the rest of this section we argue that Fp{t) < F(t), i.e., that the fidelity product is smaller than the fidelity. 
Obviously we need only diagonal elements of matrix PP^ to construct either F(t) or Fp(t), which wc write as 

N 

(pp% = |^(*)l 2 + E i p ^(*)i 2 - ( 15 ) 

For the first term we can write |Pii(t)| 2 = 1 — Oi(t) where cti(t) is some function of time with properties o»(0) = 
and < ai(t) < 1 due to relation ©. For the second term we write ^ZfcLi k^i \Fk(t)\ 2 = fiiif) where /3i(0) = and 
Pi(t) > 0. It follows that (PP^)u = l — (ai(t)—/3i(t)). By applying the trace-log formula in the same manner as before, 
we get for the fidelity F(t) = cxp(— J2i=i a i{t) + J2i=i The fidelity product corresponds to exp(— J2i=i 

which yields F(t) = Fp(t) exp(J^j = i since Xh=iA(*) > we have F(t) > Fp(t). Averaging over noise does 

not change this relation. 



IV. FIDELITY IN THE MEAN FIELD REGIME 



In this section we consider Loschmidt echo in the mean field regime, that is, by employing the Gross-Pitacvskii 
theory. The dynamics of Bosc-Einstcn condensates (BECs) is within the framework of this theory described by using 
the nonlinear Schrodingcr equation (NLSE), which wc write in dimcnsionlcss form: 

$(a>, t) + g 1D N\${x, t)\H{x, t), (16) 

where gm = ^rnX^gio /h 2 is the dimensionless coupling strength and J t)\ 2 dx = 1; here we choose giu = 0.04, 
which can be experimentally obtained by tuning the transverse confinement frequency to lj±/2tt sa 240 Hz, and with 
TV = 50. With those parameters the system is in the mean field regime with 7 rj 0.01. All parameters are identical as 
in the simulations of a TG gas, except that now uj±_ is smaller. 

To compute the fidelity of interacting BECs we repeat the same procedure as for the TG gas: first we prepare the 
condensate in the ground state of the container like potential Vl(x) (i.e. we solve numerically the stationary NLSE), 
second we suddenly expand the container to V2l(x), and solve numerically the time-dependent NLSE in the expanded 
potential without noise [V2i(x)], and with noise [V^' L (x)], with identical initial conditions. This gives us <&(x,t) and 
$'(x,£) from which we calculate the fidelity 



<9$(x,f) 



dx 2 



F GP (t) = I / $'*(x,*)$(x,t)dx| 2 . 



(17) 
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However, note that since we investigate the fidelity of a gas with N particles, the mean-field TV-particle wavefunction 
is a product state, tpGp(xi, ■ ■ ■ ,%N,t) = YljLi ®( x j: t)> an( i therefore the TV-particle mean- field fidelity is 

Fap(t) = | J i^Gpdx, . ..dx N \ 2 = [F GP (t)] N . (18) 

Finally, we average over 50 different realizations of the potential to obtain (FGp(t)) n oise and (FQ P (t)) noise- 
In Fig.^a) we plot (FGp(t)) no ise and its standard deviation for noninteracting and weakly-interacting BECs. We 
see that oscillations are superimposed on the overall decay in contrast to the TG gas case. We find that in the mean- 
field regime described by the Gross-Pitacvskii equation the fidelity decays faster for larger nonlinearity (interaction 
strength). It is worthy to point out that Fcp{t) is very dependent on the particular realization of V e (x), which is 
not the case for the TG gas. This is illustrated in Fig. [2Kb) where we show dynamics of FQ P (t) for two different 
realizations of the noise potential; we observe a large dependence of Fop(t) on a particular realization of the noise. 
This is a consequence of the fact that the oscillation frequency of fidelity |P n (t)| 2 = J2 n m l a nl 2 l a ml 2 cos[(w n — u> m )t] 
for the noninteracting BEC essentially depends on the difference between only several frequencies which is very noise 
sensitive, and this behavior is inherited in the nonlinear mean-field regime. 

Finally, in Fig. [4] we compare the fidelities of the noninteracting BEC, the weakly-interacting BEC, and the TG 
gas. Note that for proper comparison one should compare FQ P (t) with F(t). We see that the mean-field fidelity 
shows richer behavior. In the regime of parameters we used, we find that (FQ P (t)) no i se decays faster than the TG 
regime fidelity in the first part of the decay dynamics, but later the mean-field regime fidelity decay slows down in 
comparison to the TG gas decay. 
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FIG. 3: (a) Noise-averaged fidelities of evolving BECs for N = 50 and e = 0.05, and standard deviations from the noise-average. 
The averaged fidelity for a noninteracting BEC is shown with the black dot-dashed line, and its standard deviation with open 
black circles. The averaged fidelity for a weakly-interacting BEC is shown with the blue solid line, and its standard deviation 
with closed blue circles, (b) Fidelities of the evolving BECs for two different realizations of the noise potential. 
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FIG. 4: Comparison of the averaged fidelities (a) and their standard deviations (b) for the TG gas ((F(t)) no i se , red dotted 
line), the weakly- interacting BEC {(FQ P (t)) no i se , solid blue line), and the noninteracting BEC (black dot-dashed line), for the 
same number of particles. The parameters are N — 50 and s = 0.05. 
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V. CONCLUSION 

In conclusion, we have explored Loschmidt echo (fidelity) in two regimes of one-dimensional interacting Bose gases: 
the strongly interacting TG regime, and the weakly-interacting mean-field regime described within the Gross-Pitacvskii 
theory The gas is initially in the ground state of a trapping potential that is suddenly broadened, and the decay 
of fidelity is studied numerically by using a small spatial noise perturbation. We find (numerically and analytically) 
that the fidelity of the TG gas decays as a Gaussian with the exponent proportional to the number of particles and 
the magnitude of the small perturbation squared (see Fig. Q] and Fig. [2]). Our results do not depend on the details of 
trapping potential; we have obtained the same behavior for a gas that is initially loaded in the ground state of the 
harmonic oscillator potential, which is subsequently suddenly broadened. Furthermore we find that Gaussian decay 
remains if we initiate the dynamics from some excited initial state or from a superposition of such states. In the 
mean-field regime the Loschmidt echo decays faster for larger intcrparticle interactions (nonlinearity) , and it shows 
richer behavior than TG Loschmidt echo dynamics with oscillations superimposed on the overall decay (see Fig. [3] 
and Fig. [4]); it also has much larger sensitivity on the noise (see Fig. [3]f b) ) . Finally, we would like to mention that 
perhaps the most interesting regime of Loschmidt echo dynamics would be for intermediate Lieb-Liniger interactions, 
which seem to be exactly solvable only for specific external potential configurations (2(| . 
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